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Abstract 

We relate the U-duality invariants characterizing two-center extremal black hole solutions in the 
stu, st^ and t^ models of = 2, d = 4 supergravity to the basic invariants used to characterize 
entanglement classes of four-qubit systems. For the elementary example of a D0D4-D2D6 composite 
in the model we illustrate how these entanglement invariants are related to some of the physical 
properties of the two-center solution. Next we show that it is possible to associate elliptic curves 
to charge configurations of two-center composites. The hyperdeterminant of the hypercube, a 
four-qubit polynomial invariant of order 24 with 2894276 terms, is featuring the j invariant of the 
elliptic curve. We present some evidence that this quantity and its straightforward generalization 
should play an important role in the physics of two-center solutions. 

PACS numbers: 03.67.-a, 03.65.Ud, 03.65.Ta, 02.40.-k 
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I. INTRODUCTION 



The aim of the present paper is to show that it is possible to relate the entanglement 
measures usually used in studies concerning four-qubit systems to the U-duality invariants 
found recently by Ferrara et. al.- characterizing extremal two-center black hole solutions in 
the stu, st^ and models. Interestingly as a byproduct of establishing this correspondence 
one can also come across an interesting connection between charge configurations of such 
black holes and a special class of elliptic curves. 

Multi center black hole solutions provide an interesting research direction within the 
rapidly evolving field of black hole solutions in supergravity, string and M-theory. For 
such solutions the attractor mechanism^ has been generalized giving rise to split attractor 
fiows'^'^. For two-center solutions the latter term refers to the situation, when in moduli 
space after crossing walls of marginal stability the attractor fiows are separately evolving to 
the attractor points of the constituent single center solutions. Recently these developments 
have triggered activity in a variety of new research fields such as attractor fiow trees, entropy 
enigmas, microstate counting, and bound state recombination^""—. 

For single center solutions it has become obvious that the notion of duality charge orbits 
and their invariants^i are useful concepts for classifying black hole solutions together with 
their supersymmetry properties. Except for a special case^^ for multicenter solutions the 
corresponding structure of orbits and invariants is still unknown. In order the remedy this 
situation in a recent paper Ferrara et.al. conductectiii^ a systematic study on the structure 
of invariants characterizing the charge configurations and invariants of two center solutions. 
In the case of the stu^ st^ and models looking at the structure of such invariants one 
immediately notices structural similarities to the well known sets of four-qubit invariants^ii^ 
discussed in the seemingly unrelated field of quantum information. 

Since the advent of the black hole qubit correspondence^^ such coincidences should not 
come as a surprise. It is well-known by now that few qubit entangled systems are capable of 
providing interesting new insight into the structure of black hole solutions and their attractor 
fiows^^. The occurrence of these qubits is related to the presence of tensor products of 
the spin ^ irreducible representations of SL{2) groups. Such products of S'L(2)s show up 
as subgroups of [/-duality groups governing the entangled web of dualities of supergravity 
models giving rise to black hole solutions. Initially mathematical coincidences were recorded 
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merely for two and three qubit systems and the corresponding axion-dilaton and stu black 
holes^, however evidence for n-qubit systems with n > 3 to make their presence in this 
context started to accumulate. 

In this line of development the relevance of four qubit systems to black hole solutions in 
supergravity was first pointed out in our paper--?- where the isomorphism 

so(4, 4) ~ s/(2)^© (2,2,2,2) (1) 

has been used to describe 7 x 16 of the 133 -£'7(7) generators of = 8, c? = 4 supergravity. 
These generators describing seven copies of four-qubit states are not belonging to the s/(2)^ 
subalgebra of -^7(7). A suitable incidence geometry accounting for the relationship between 
the 7 groups of 16 generators is that of the dual Fano plane giving rise to a geometry, dual 
to the one describing the "tripartite entanglement of seven qubits" interpretation2Si2i of the 
quartic -^7(7) black hole entropy formula. The isomorphism of Eq. ([T]) was also discussed in 
the review paper of Borsten et. al.— providing further interesting examples of simple qubit 
systems. 

In our next paper>^ in a four-qubit entanglement based formalism the structure of ex- 
tremal stationary spherically symmetric black- hole solutions in the STU model of = 2, 
d = 4 supergravity was described. The basic idea facilitating this interpretation was the fact 
that stationary solutions in d = 4 supergravity can be described by dimensional reduction 
along the time direction^i. In this d = 3 picture the global symmetry group SL{2, M)^ of 
the STU model is extended by the Ehlers SL{2, M) accounting for the fourth qubit. One 
can then introduce a four-qubit state depending on the charges, the moduli, and the warp 
factor. Here it was also noticed that in the terminology of four-qubit entanglement extremal 
black- hole solutions should correspond to nilpotent, and nonextremal ones to semisimple 
states. The upshot of these considerations was the emerging possibility of relating the en- 
tanglement properties of such and similar states to different classes of black-hole solutions 
in the STU model. The challenge of elaborating on this idea was recently taken up in the 
papers of Borsten et.al.— In these papers the authors applied the black-hole qubit correspon- 
dence to the problem of classifying four-qubit entanglement. The key technical ingredient 
was the Kostant-Sekiguchi theorem which establishes the link between nilpotent orbits of 
extremal black holes and four qubit entanglement types. The emerging picture is: we have 
31 entanglement families which reduce to nine up to permutations of the qubits. These nice 
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papers confirmed once again that the input coming from string theory can be useful in estab- 
lishing results in a different field, since the literature until now on four qubit entanglement 
classification was confusing and seemingly contradictory. 

In this paper we would like to show that the charge orbit classification of two-center black 
hole solutions in the stu model is another arena where four-qubit systems naturally make 
their appearance. As a first possible step in this direction here we establish a correspondence 
between the U-duality invariants of Ferrara et.al.- and the four-qubit invariants showing up 
in classification schemes of entanglement types in quantum information. Establishing this 
correspondence simplifies some of the invariants proposed so far, clarifies their geometric 
and algebraic roles, and provides hints for further generalizations outside the framework 
of N = 2, d = 4 supergravity. As an extra bonus the four-qubit picture also hints at a 
basic physical role these invariants are playing in the theory of two-center solutions. For 
one of the invariants not fully appreciated yet, our considerations establish a special role. 
It is the SL(2)^^ and permutation invariant hyperdeterminant of type 2x2x2x2. This 
is a polynomial of order 24 in the 16 amplitudes of the four-qubit state. Mapping the 
16 amplitudes to the 16 charges characterizing two-center solutions in the stu model, for 
a special case we show that the structure of this hyperdeterminant seem to govern issues 
of consistency in the realm of two-center solutions. These ideas also suggest a natural 
way for associating an elliptic curve of a special kind to a particular charge configuration. 
The coefficients of our elliptic curve are the algebraically independent four-qubit invariants, 
and its discriminant is just the hyperdeterminant. We also present some evidence for the 
conjecture that the structure of the j invariant of the elliptic curve should play an important 
role in the physical properties of the two-center solution. The idea that elliptic functions 
and the j invariant might play some role in four-qubit systems and the black hole qubit 
correspondence was first suggested by P. Gibbs^^ some related discussion appeared in the 
paper of Bellucci et.al.— . 

The organization of this paper is as follows. In Section II. we summarize the background 
material on four-qubit invariants, reduced density matrices, and the structure of the hyper- 
determinant of the hypercube. We introduce a quartic polynomial featuring the algebraically 
independent four-qubit SL{2, C) invariants. In Section III. we are discussing extremal two 
center black hole charge states in a four-qubit based picture. Here we work out a dictio- 
nary between the invariants found by Ferrara et. ah" in the so called Calabi-Visentini basis 
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and the algebraically independent four-qubit ones in the "special coordinates" basis. We 
give some of the invariants used in Ref.- a simpler appearance, and connect other invari- 
ants of physical meaning to properties of four-qubit reduced density matrices. Here we also 
show that the set of algebraically independent polynomial invariants in both the four-qubit 
and the Ferrara et.al. description are based on two seemingly different quartic polynomials 
however, with the same resolvent cubic. 

Section IV. is devoted to a case study featuring BPS DODA — D2D6 composites in the 
model. In the paper of Bates and Denef^ this elementary example has already turned out 
to be a good playing ground for investigating the basic properties of two center solutions, 
hence we opted for illustrating the physical role of our four qubit invariants in the very same 
setting. These considerations relate the (necessary) consistency condition, guaranteeing the 
BPS composite to exist, to the positivity of the hyperdeterminant and to the extra constraint 
that the two nonzero invariants of the model are having the same sign. It turns out that 
precisely these conditions are the ones guaranteeing the fundamental quartic polynomial to 
have real roots. We then associate an elliptic curve of Weierstrass canonical form to the 
resolvent cubic of this quartic and show how physical properties are nicely encapsulated in 
the structure of its j invariant. 

In Section V. we examine the status of our rather ad hoc assignment: two-center charge 
configuration-elliptic curve, more thoroughly. By switching to the most general Tate form of 
an elliptic curve we show that our association of elliptic curves to two-center black hole charge 
configurations in the stu model is a natural one. This means that the nonzero coefficients 
aj with j = 1,2,3,4 appearing in the Tate form are algebraically independent four-qubit 
invariant homogeneous polynomials of order 2j. These coefficients have important physical 
meaning: ai is just the canonical symplectic pairing between the charge vectors of the two 
centers. Vanishing of gives rise to the st"^, and a further vanishing of 02 results in the 
truncation. One can also see that in the st^ and models sending 03 to zero corresponds 
to the limit when our elliptic curve degenerates. One of the nontrivial coefficients in the 
Tate form is always zero for the stu model. This is related to the vanishing of a nontrivial 
polynomial constraint of homogeneous degree 12 valid in the stu model, already observed 
by Ferrara et.al.i Based on the latest results of Andrianopoli et. alM we conjecture that we 
should be able to generalize our correspondence between charge orbits and elliptic curves 
also for the case of maximal N = 8, d = 4 supergravity. In this case the stu model should 
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arise as an ae = truncation implemented by the vanishing of a polynomial of order 12. 

The aim of our last speculative Section VI. is to draw the readers attention to some 
interesting structural similarities showing up in a variety of physical contexts where our 
four-qubit invariants parametrizing elliptic curves might play a crucial role. Here we give 
a new look and interpretation to a triality invariant curve originally introduced by Seiberg 
and Witteni^. Now this curve is paramaterized by four-qubit invariants also displaying 
permutation invariance. In this new setting we also invoke the F-theory interpretation of 
this curve as was given by Ser>i^. Finally our conclusions and some comments are left for 
Section VI. 

II. FOUR QUBIT SYSTEMS 

In order to facilitate a four-qubit description of the two center charge configurations in 
the stu model our aim in this subsection is to review the background material on four-qubit 
states and their entanglement measures. A four qubit state can be written in the form 

I A) = 5Z ^iohi2i3 No«1^2«3) , \ioHi2i3) = Ko) ® Kl) ® K2) ® Ks) e K) (8) 14 ® ^2 ® V3 (2) 

where Vo,i,2,3 = C^. Let the subgroup of stochastic local operations and classical 
communication^^ representing admissible fourpartite manipulations on the qubits be just 
SL{2, C)®^ acting on |A) as 

\A)^iSo®Si0S2 0S3)\A), S^eSL{2X), a = 0,1, 2, 3. (3) 

Our aim is to give a unified description of four-qubit states taken together with their SLOCC 
transformations and their associated invariants. As we will see states and transformations 
taken together can be described in a unified manner using the group 5*0(8, C). 

Let us discuss the structure of four-qubit SL{2,C)^'^ invariants^'i^'^i^. The number of 
algebraically independent four-qubit invariants is four. We have one quadratic, two quartic, 
and one sextic invariant. In our recent paper— we investigated the structure of these invari- 
ants in the special frame where two of our qubits played a distinguished role. As we will see 
this scenario is just the one needed in the two-center STU black hole context since in this 
setting one of the special qubits (the one labelled by the number 0) will be associated to the 
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horizontal SLh{2,M.) of Ferrara et.al.i and the other (the one labelled by the number 1) is 
arising as the first factor from the structure SL{2, M) x 5*0(2, 2) known from the STU model. 
Indeed such structure is the one arising as a special case of the infinite Jordan symmetric 
sequence of = 2 = 4 supergravity theories^i. 

To an arbitrary state |A) we can also associate a 4 x 4 matrix 

/Annnn Annm Anmn AnmA f A^ A A 

^1 ^2 ^3 ^4 

Ci 



c 



^0000 Aoooi AooiO Aooil^ 

^0100 AqioI Aquo Aqih 

Alooo Alool Alolo Aioii 



\Aiioo Aiioi Aiiio Aiiiiy 



(4) 



or four four- vectors. The splitting of the amplitudes of |A) into such four- vectors reflects 
our special choice for the distinguished qubits compatible with our conventions. We will 
also need the matrices 



M 



(a^ A^ B' 

(Jl Q2 jjl Jj2 
^3 ^4 ^3 ^4 



(a^ A^ B^\ 
A^ A^ B^ B^ 

\C2 d2 d^J 



(5) 



Notice that in four-qubit notation the matrices Ai and Af are arising from C by the permu- 
tations (012) (3) and (0)(123) meaning that the index structure of these matrices is 

Notice that these matrices are entering in the reduced density matrices as 

^oi = Tr23|A)(A|=££t = Tro3|A)(A| = TWA^^ ^02 = Tri3|A)(A| = A^A^t (7) 

^23 = Troi|A)(A| = /:t£ ^03 = Tri2|A)(A|=A^t_^ ^13 = TrosIA) (A| = A/'+A/". (8) 

where overline denotes complex conjugation. 

Now we introduce on the vector space C'^ x corresponding to the third and 

fourth qubit a symmetric bilinear form g : x — )■ C with matrix representation 



g = 6 ®e 



l\ 




^0 1\ 


[-1 0) 




i-i 0) 



(9) 
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This means that we have an SL{2, C)^^ invariant quantity with the exphcit form 
g{A, B)=A-B = QapA^'B'^ = A^B"" = A^B^ - A^B^ - A^B"^ + A^B\ 
We can also introduce a dual four-qubit state 

|A) 

with the associated matrix 



(10) 



«0«l*2i3=0,l 



(11) 







■^0001 


Aqoio 


Aooii^ 


/a' 






/ = 


-^0100 


-^0101 


Aoiio 


Aoiii 


b' 


b^ 


b^ b" 


A 1000 


Alool 


Alolo 


Aioii 










V^llOO 


AllOl 


Aiiio 


Aiiii/ 






d^j 



(12) 



where 

a» = e''f^''^BpC^Ds, b^ = e^^^^^A^C^Ds c'' = e""^^^ A^BpDs d^ = e"^^^ A^BpC^. (13) 

Here e^^'^^ = +1, and indices are lowered by the matrix of g. Notice that the amplitudes of 
the dual four-qubit state are cubic in the original ones. 

Using these definitions we define the quadratic and sextic invariants as 



h = \{A-D-B-C), h = h^a-d-b-c). 



(14) 



The explicit form of the sextic invariant in terms of the dot product of Eq.f llOp is 

^ A- A A- B A-d\ / a . n n . n n .r<\ 

A-C B -C C ■ D -Det 
A-D B-D D-D 



2h = Det 
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A-B B-B B-C 
A-C B-C C-C 
A-D B-D C-D 



(15) 



We also recall that the explicit form of Ji is hiding its permutation invariance. Moreover, 
though the expression of I3 of Eq.( fT4l) is similar to the one of Ii the invariant I3 is not 
invariant under the permutation of the qubits. 

Now we turn to the structure of quartic invariants. We have two independent of such 
invariants^ and the simplest of them is the obvious expression 



Det£ 



(16) 
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i.e. the determinant of the 4x4 matrix of Eq.(|l]). In order to present the definition of the 
second one we define separable bivectors of the form 

li^uaP = ^f,aKl3 - ^^ipKa, «, /3, /i, Z/ = 1, 2, 3, 4. (17) 

Here our labelhng convention A^q, indicates that /i = 1,2,3,4 identifies the four-vector in 
question (i.e. A^B^C or D of Eq.(jl])), and the label a = 1,2,3,4 refers to the component 
of the particular vector. Now our last invariant is the quartic combination 

h = ^U^ua^W^"''''. (18) 

Obviously the symmetric nondegenerate bilinear form g of Eq.f llOl) acting on four vectors 
like A, B,C,D G is inducing a corresponding symmetric nondegenerate bilinear form on 
the space of bivectors /\^C^. By an abuse of notation we use again the symbol ■ for this new 
bilinear form with the definitions^ 

{AAB)-{CAD) = 2{{A ■ C){B ■ C) - {A ■ D){B ■ C)). (19) 

Now I2 can also be written in the equivalent form 

1 .1 J 

(20) 



^ 6 



{A A B) ■ {C A D) + {A AC) ■ {B A D) - ]^{A A Df - ^(5 A Cf 



An important comment here is in order. Let us have a look at 1^ and also at the deter- 
minants of the matrices of Eq. ([5]) 

L = h = Det£, M = DetM, N = DetAf. (21) 

Then one can prove^*^ 



One also has the constraint 



L + M + N = 0. (22) 
M - AT = 3/2 - 2/2, (23) 



that we will need later. 

It is knowni^ that the minimal set of algebraically independent SL(2)^'^ invariants is 
consisting of a quadratic, two quartic and one sextic invariant. Our choice for this set will 
bei^: Ii, 12,14 and J3. Let us now present the reason for this choice. Let us consider the 
matrix 

n = CgC^g. (24) 



Then its characteristic polynomial is 

T.^{Ki,i,^,,t) = Det(lt ~n)=t^- + Qht" - 4/3t + ll (25) 
Clearly by Newton's identities we have 

/i = ^Trfi, /2 = ^[(Trl])2-™2], (26) 

Jg = ^[(Trn)^' - + 2™^], (14)^ = Detfi. (27) 

This form of writing our invariants is related to the fact that there is a 1 — 1 correspondence 
between the SL{2, C)^^ orbits of four-qubit states and the 50(4, C) x 50(4, C) ones of 4 x 4 
matrices. 

The polynomial of Eq. f|25p in the four-qubit context appeared in our recent paper^- its 
role as a characteristic polynomial has been emphasized in Ref.— . The discriminant of this 
fourth order polynomial is the hyper determinant^ D4 of the 2x2x2x2 hypercube AjQj^jjjg. 
It is a polynomial of degree 24 in the 16 amplitudes and has 2894276 terms^. D4 can be 
expressed^ in terms of our fundamental invariants as 

256^4 = - 27T^ (28) 

where 

S = {II - II) + A{ll - hh), T = {II - - h) + (/3 - hhf. (29) 

For an alternative form of D4 see the papers of Refs.— i^^. 

In closing this section we briefly discuss some results on the classification of entanglement 
classes for four qubits^i^i. By entanglement classes we mean orbits under SL{2, C)^^ ■ Sym4 
where Sym4 is the symmetric group on four symbols. The basic result states that four qubits 
can be entangled in nine different waya^"^. It is to be contrasted with the two entanglement 
classes^^ obtained for three qubits. For a refined classification of four qubit entanglement 
motivated by the black hole qubit correspondence see the papers of Borsten et.al.— 

Let us consider the matrix 

= f ° (30) 
\-A^g J 

which now can be regarded as an element of the Lie algebra of 50(8, C). If the matrix TZ\ 
is diagonalizable under the action 

1 I 5o ® 5i 

7^A^57^A5-^ 5= , 5„g5L(2,C) 31 

52® 
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we say that the corresponding four-qubit state |A) is semisimple. If TZa is nilpotent then we 
call the corresponding state |A) nilpotent too. It is known that a nilpotent orbit is conical 
i.e. if |A) is an element of the orbit then A|A) is also an element for all nonzero complex 
numbers A. Hence a nilpotent orbit is also a G'L(2,C)^^ orbit. It is clear that for nilpotent 
states all of our algebraically independent invariants are zero. 

A semisimple state of four qubits can always be transformed to the forn>2i 

\Gatcd) = ^(10000) + 11111)) + ^(10011) + 11100)) 

+ ^(10101) + 11010)) + ^(10110) + 11001)), (32) 

where a, b, c, d are complex numbers. This class corresponds to the so called GHZ class 
found in the three-qubit case2^. For this state the reduced density matrices obtained by 
tracing out all but one of the qubits are proportional to the identity. This is the state with 
maximal four-partite entanglement. Another interesting property of this state is that it does 
not contain true three-partite entanglement. A straightforward calculation shows that the 
values of our invariants (Ji, I2, 13, 14) occurring for the state \Gabcd) representing the generic 
class are 

I, = ha^ + b^ + c^ + d% h = l[{ahf + {acf + {adf + {hcf + {hdf + [cdf], (33) 
4 o 

Ir^ = ^[{ahcf + {ahdf + {acdf + {hcd)\ h = abed (34) 

hence the values of the invariants (4Ji, 6J2, 4J3, J|) are given in terms of the elementary 
symmetric polynomials in the variables (ti, ^2, ^3, ^4) = (a^, 6^, c^, c?^). For the semisimple 
states I Gabcd) the value of can be expressed as^'^ 

A = ^n,<,(t,-t,)2, {ti,t2,h,U) = {a\b^,c^d^). (35) 

Notice that for the states \Gabcd) with D4 nonvanishing (tj 7^ tj) the corresponding matrix 
of Eq. (!30|) belongs to a Cartan subalgebra of 5*0(8, C). The stabilizer of such states 
corresponds to the Weyl group of SO {8, C). This stabilizer is the Klein group generated by 
the four elements / / (g) / (g) / and aa <^ Ca ^ (Ja ® era for a = 1, 2, 3. 

III. TWO CENTER EXTREMAL BLACK HOLES AS FOUR QUBITS SYSTEMS. 

In the paper- of Ferrara et.al. in order to describe the structure of the U-duality invariant 
polynomials associated to the two-center extremal black hole solutions the Calabi-Visentini 
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(CV) basis has been used. Here by [/-duality we mean the continuous hmit- vahd for 
large values of the charges of the usual-nonperturbative string theory symmetries. First we 
describe the connection of the CV basis to the one making the four-qubit structures explicit. 
Next we turn to an entanglement based understanding of the structure of the two-center 
[/-duality invariants. 

In the CV basis the two-center black hole solutions are characterized by a pair of real 
charge vectors Qi, Q2 G 

Qi = {P^,P\P^,P\Qo,QuQ2,Q3f, Q2 = {p'',p\p',p',qo,qi,q2,q3f. (36) 

As we see these charge vectors are containing two four- vectors each namely P^ , Qi = 
rjijQ'^ , and p^,qi = rjuq-^ I,J = 0, 1,2,3 where the raising and lowering of the indices / 
and J are effected by the metric rjjj and rj^'^ of 5*0(2,2) answering the symmetric bilinear 
form h acting on the charge four-vectors P^ and Q"' in the CV basis as 

h{P, Q)=PoQ = i^jjP'Q' = -P'Q' - P^Q' + P^Q^ + P^Q\ (37) 

Let us now relate the 16 component charge vector in the CV basis characterizing a 
particular two center extremal black hole solution in the STU model to a real unnormalized 
four-qubit pure state |A) by relating the four-vectors P^,Q^,p^,q^ G of Eq. (!36l) and 





of Eq.g 


1 as follows. 












/ - \ 












1 


A^ + A^ 




Q' 


1 


B^ + B^ 


p2 




-A^ - A^ 






"7! 


-B^ - B^ 


\P' J 


CV 


\-A^ + 






CV 


\-p2 + By 


fp^ \ 












^ D'-D^\ 


pi 


1 






q' 


1 




P' 


"71 






q' 


"7! 


-D^ - 


\P' J 


CV 


+ C7 






CV 





(38) 



(39) 



Then using the definitions above we clearly have for example 

PoQ = A-B = g{A, B) = Qo^pA'^B^ (40) 
with the bilinear form g defined as in Eq. (llOI) . 
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Let us also give the connection between the Calabi-Visentini basis and the one usually 
used in special geometry, i.e. the special coordinates (SC) symplectic frame. This frame 
yields the usual set of electric and magnetic charges i.e. {P^,Qi)sc and {p^,qi)sc- In the 
following we will use these charges so it is important to clarify their relationship to the 
components of our four qubit state |A) (see Eq.(jl])). 





















( ^ 


^0001 








p2 




Aoioi 








Aoolo 




A' 




p3 




A-ollO 






Q2 


VAooii/ 




Uv 






SC 


VAoni/ 




Uv 





(41) 















( Alloo^ 








Aiooi 












Aiioi 






gs 


Aioio 








pi 




Aiiio 






(I2 


VAioii/ 










SC 


VAiiii/ 




Uv 


\-1o) 



(42) 



We note here however, that our conventions are slightly different from the ones used in Ref.l. 
The charges p^,p'^,p^, go and P^, P^, P^, Qo in the SC basis used by us are the negatives of 
the corresponding ones in the SC basis as used in Ref.l. see Eqs. fl38ti39|) . Eqs. (1411142]) and 
Eq. (B.3) of that paper. 

Notice also that in our four-qubit state |A) with amplitudes Ajgj^jjjg sit two three-qubit 
states with amplitudes Aoiiiais! and Aijjjjjg. The first set of amplitudes labelled by zq = 
corresponds to the charge configuration of the first black hole and the second labelled by 
zo = 1 describes the second black hole. As we see the first label plays a distinguished role 
with an extra SL{2, M)o (dubbed by Ferrara et.al. the horizontal one) acting on. This group 
represents the generalized exchange symmetry between the centers. 

Let us now define three bivectors 

X = AAB, Y = CAD, Z =^{AAD - B AC). (43) 

In component notation we have for example 

Xap = A^Bp - A^Ba. (44) 

Switching to the Calabi-Visentini basis these objects are the T-tensors used in the paper of 
Ferrara et.al.- As we can see the bivectors X and Y are separable i.e. they are precisely the 
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ones satisfying the Pliicker relations 

-^12-^34 — -^13-^23 + -^14-^23 = 0, ^12^34 " ^13^23 + ^14^23 = (45) 

on the other hand the bivector Z is entangled i.e. in the nomenclature of fermionic 
entanglement^- it has Slater rank two. 

Now we introduce the shorthand notation for the product of two separable bivectors as 
defined in Eq.([l9l) 

X -Y = {A^B)■{C ^D) (46) 
Notice that in the notation of Eq. (HT]) -f H2]) 

= X-X = 2{A^B^ -{A-Bf), = Y -Y = 2{C^D^ - {C ■ Df) (47) 

are just two times the quartic invariants of the charges characterizing the two black holes 

\x^ = hiQi) = -D■s{Aon^,^,), \y^ = hm = -/^3(Ain...3)- (48) 

Here D^^Ai-^i^i,^) is Cayley's hyperdeterminant^. 
With these definitions we can define the quantities 

/+2 = ^X^ /+i = ix-Z, Io = ^{2Z'-X-Y), U=^-Y-Z, I-2 = \y^ (49) 
and the ones 

I' = \x- Y, I" = ^Z\ (50) 

When reinterpreted in the CV basis these are precisely the SL{2, R)i x SL{2, R)2 x SL{2, R)^ 
invariants of Ferrara et.al.- The important property of the invariants of Eq.( H9|) is that they 
are covariants with respect to SL{2, R)q acting on the distinguished (horizontal) qubit. 
Indeed, they are sitting in the 5 (spin 2) irreducible representation of this group. 

Now we elucidate another aspect of this important property of the set of invariants of 
Eq.f HOj) . In order to do this we first recall that for two-qubits the canonical measure of pure 
state entanglement is the concurrence^ 

C = 2\D2{A^,^,)\ = 2|AooAn - AoiAiol (51) 

which is related to the determinant of an ordinary 2x2 matrix. For three-qubits the basic 
quantity characterizing genuine three-qubit entanglement^ is the three-tangle 

T = A\D,{A,,,,i,)\ (52) 
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where now D-^ is Cayley's hyperdeterminant^i^^. According to the method of Schlafli^^ 
is related to the discriminant A2 of the quadratic polynomial 

^2(Aio^l^2, ^) = D2i^0^^^2t + Al^^J = (Aq " Ao)t' + 2(Ao " Al)t + (Ai ■ Ai) (53) 

where Aq and Ai are four-vectors with components (Aqoo, Aqoi, Aqio, Aqu) and 
(Aioo, Aioi, Alio, Am) and the ■ product is the usual one of Eq. fnU]) . Moreover, due to 
permutation invariance of we obtain the same expression whenever the first or the sec- 
ond qubit plays a distinguished role. Obviously the quantities J+i = (Aq-Aq), Jq = (Aq-Ai), 
and J_i = (Ai • Ai) are 5*27(2)1 x SL{2)2 invariants, however they are covariants under the 
"horizontal" SL{2)o. Indeed, the triple J+i, Jq, J_i transforms according to the irreducible 
representation 3 of spin 1 of this "horizontal" group which can be regarded as some sort of 
generalized exchange symmetry working between the two two-qubit systems. It is also clear 
that searching singlets with respect to this "horizontal" symmetry group can reveal some 
new properties of our pair of two qubit systems, namely that they are secretly comprising a 
system having a higher degree of symmetry. 

Now proceeding by analogy we define the polynomial 

n4(A,oi,i,,3, t) = D3{Ao^,^,i,t + Au,i,i,) = 1+2^^ + 4J+it=^ + Qlot^ + 4/_it + /_2. (54) 

A straightforward calculation shows that the coefficients of this polynomial are precisely 
the covariants of Eq. (H9|) . Now according to theorem 14.4.1 and corollary 14.2.10 of Ref.— 
the discriminant of this quartic polynomial A4 divided by 256 is just the hyperdeterminant 
DA{Aigi-^i,ji.^) of the hypercube of format 2x2x2x2. Besides being a singlet with respect 
to the horizontal SL{2)q, D4 is also an invariant under 6*4 the permutation group of the four 
qubit system. 

At this point one can notice-^ that the polynomials S4 and 114 of Eqs. fl25p and f l54p are 
having the same discriminants A4 hence both can be used to obtain an expression for D^. 
Notice that S4 is a polynomial with its coefficients Ii, I2, 13, 14 also being SL(2)q singlets 
however, 114 is a polynomial with coefficients I±2,I±i,Io being merely SL{2)o covariants. 
Using this observation we can construct new SL{2)q singlets from the quantities of Eq. (H^ 
by relating them to the known algebraically independent four-qubit ones namely Ii, 12,13 
and I 4. 

In order to relate the SL{2)o singlets found by Ferrara et.al. to our four algebraically 
independent four-qubit invariants one just has to compare the relevant expressions. In fact 
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many of their invariants and the constraints satisfied by them take in this four-qubit setting 
a much simpler and instructive form. In particular their complete set of invariants with 
corresponding degrees 2, 4, 6, 8 denoted by W, x, Iq, Tr(j7'^) is related to ours as 

W = 2Ji, x = 3/2 -2/2, h = -l3- (55) 

By virtue of Eqs. (123|) and (150|) we also have the relations 

I'-I" = ^l2, x = M-N. (56) 

We still have to account for the invariant Tt{J''^) of order 8 built from J' the symmetric 
traceless matrix comprising the covariants /±2, -^±1; -^o- The 5 independent components of 
this matrix are transforming according to the 5 of 5*17(2)0. In order to reveal the meaning 
of this invariant and also an extra one Tt{J'^) of order 12 let us reconsider the awkward 
looking polynomial constraint of Eq. (5.6) of Ref.l. 

V,, . P^^W^I.^M^'H^ + + + JL . 0. (57) 

Using the dictionary of Eq. (l55l) we can cast this constraint in the nice form 

Tr(j3) = [4(/2 _ 1^/3) _ Tr(j2)](/2 _ j^) _ (J3 _ j^j^)2_ (53) 

There is one more invariant^ of order 8 which is directly related to Tt{J'^) 

Vs = -12Tr( + 24/6>V + (W' + 2xf. (59) 

Recall that the constraint Vs = implements the reduction of the stu model to the st'^ 
modeU- in a manifestly SL{2)q invariant manner. Using again Eq. fl55|) the new form of Vs is 

Vs = 12(3/2 - 4/1/3 - Tt{J^)). (60) 

Putting this into Eq. (l58|) and recalling Eq.f l29|) one obtains the simple expressions 

Tt{J^) = S, TiiJ^) = -T, (61) 

provided 

Vs = = -12L^. (62) 
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The first result of tliese considerations is tliat — j^'Ps is really the square of the basic 
fourth order invariant L = 14^. According to Eqs.(|l]) and fl2T]) L is just the determinant of 
the matrix C we have started our four-qubit considerations with. Moreover, according to 
Eqs.(I7HH]) we also see that for unnormalized four-quhit states we have 

- Y^Vs = Det^oi = Det^i2, (63) 

hence this invariant is related to the determinant of the reduced density matrices of our 
four-qubit state |A) corresponding to the bipartite split of the form: (01)(23). We can also 
conclude that in the four-qubit picture the reduction of the stu model to the st"^ is effected by 
sending one of the eigenvalues of the reduced density matrices corresponding to the (01) (23) 
split to zero. Notice however, that the remaining density matrices corresponding to the 
remaining two splits (02) (13) and (03) (12) are generally not sharing this property. This 
means that for the st"^ model we have 

L = 0, M^O, N^O. (64) 

Recall now that a suitable further reduction to the model is obtained by employing 
the following two SL{2)q invariant constrainta^ 

X = 0, Vs = 0. (65) 

By virtue of Eq. ( l56l) and the constraint L + M + = these constraints can be described 
in the compact form 

L = M = N = 0. (66) 

This means that for the reduction of the stu model all of the reduced density matrices of 
Eqs.([7}l8]) of the four qubit state |A) have a zero eigenvalue. 

The second result of our considerations finally clarifies the role played by the invariants 
Tr{J^) and Tt{J^). In particular looking at Eqs.(l2HD, (EH]) and Eq.(^ we see that that the 
four-qubit hyperdeterminant 256D4 which is just the discriminant of our polynomial S4 of 
Eq.( l25|) can be expressed with the help of these invariants of order 8 and 12 as — 27T^. 
Moreover, since the discriminant of a quartic is the same as the discriminant of its resolvent 
cubic one can also show that the two quartic equations S4 = and 114 = of Eqs. and 
(IMll are having the same resolvent cubics. Indeed, a straightforward calculation shows that 
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the corresponding resolvent qubics in both cases are of the form 

u^-Su-2T = 0. (67) 

In the first case we get back to the known expressions for S and T of Eq. (!29|) . and in the 
second one we get^ 

S = Ml - AI+J^i + /+2/-2, T = I^ + + /'1/+2 - /+2/0/-2 - 2/+1/0/-1. (68) 

IV. THE D2DG - 0004. SPLIT 
A. Invariants 

In order to uncover the role of our four-qubit invariants playing in the theory of two 
center black hole solutions let us consider a special class of two center black hole solutions 
in the model featuring a DODA — D2DQ split in the type HA duality frame. In this case 
the vectors of Eq. ( l4Tl) - (l42l) in the special coordinate (SC) basis are 

_ P 
~ P ' 

corresponding to the first black hole with charge configuration Qi of a BPS DODA system 
and 

corresponding to the second black hole with charge configuration Q2 of a BPS D2D6 system. 
For BPS configurations in the first case we should have 

- DsiQi) = {A ■ A){B ■B)-{A- Bf = AUP^ > (71) 
and in the second the corresponding constraint is 

- ^3(22) = {C ■ C){D ■D)-{C- Df = -Avq" > 0. (72) 
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Notice that our charge spht for the special values of?7 = 4, P = g = l and v = —4 
incorporates the illustrative example of Bates and Denef^ (in that paper v is related to ours 
via a sign flip). The four- vectors A", B^, and are comprising the 16 amplitudes of a 
real four-qubit state as displayed in Eqs. and Now using Eqs. ffT^ - ffT^ . ffTB]) and 
(12Uj) the algebraically independent four-qubit invariants Ji, J2,/3 and J4 can be calculated. 
The explicit forms of these invariants are 

h = ^{Uv-3Pq), l2 = lll h = -Pq{Pq + Uv)\ h = 0. (73) 

Notice that by virtue of Eqs. fl22|) - fl23|) our DODA — D2DQ example illustrates the con- 
straints we have already discussed in connection with the model, namely the ones 
L = M = N = 0. Since L = I4 = and 3/2 = 2If {x = 0) the quartic equation 
S4 = arising from the polynomial of Eq. (J25l) now reduces to a cubic one of the following 
form 

- Aht^ + Allt -4/3 = (74) 
The discriminant of this cubic equation is 




4/3 I /3-^^) |. (75) 

According to Eqs. fl28p -f l29p the hyperdeterminant D^^ is also related to this discriminant 
and is of the form 




256D4 = 27/| \[^] - /3 = -{UvPq){PqY{Uv + 9Pqy{Uv + Pqf. (76) 



B. Consistency condition 

Let us now consider the necessary condition's^ for our two center charge configuration 
supporting a corresponding two center stationary extremal BPS black hole solution. As it 
is well-known this condition is of the form 

|x,-x2| = (Q„Q2) Jf (77) 
2Im(Z2Zi),.=oo 

where xi 2 are the locations of the centers, r = |x|, Z12 are the central charges corresponding 
to the charges Qi,2, and the symplectic product of the charge vectors (Qi, Q2) is related to 
our quadratic four-qubit invariant as 

/i = ^(Qi,Q2). (78) 
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The explicit forms of the central charges for our centers are 

Here r is as usual the complex scalar field of the model 

T = X - iy, y > (80) 

with negative imaginary part^"^, and K = — log(8|/^) is the Kahler potential. 
Now the constraint dictated by Eg. (1771) is 

/iIm(Z2Zi)oo > 0. (81) 

Explicitly we have 

8yllm(Z,Z,)^ = -3Pv{xl + ylf + (9Pg + Wv)xl^ + (9Pg - Uv)yl + Wq. (82) 

We are interested in the structure of BPS DODA — D2D6 composites hence according to 
Eqs.(I7I])-(I72]) we should have 

UP > 0, vq< 0. (83) 
For all possible sign combinations satisfying these constraints we have 

{-Pv)h < 0. (84) 
Let us now write our consistency condition as 

(-k) (-^-I-(^2^i)oo) - (-^) ^ > 0. (85) 
Now a calculation shows that by virtue of Eg. (1841) consistency demands that 

V = 3 {Pv{xl, + yl) + h-UvY + UP{2vy^f -\{Uv + Pq) {Uv + 9Pg) < 0. (86) 

Let us now look at the expression for our hyperdeterminant of format 2x2x2x2 as 
given by Eq. fl76p . Clearly positivity of implies that —UvPq > and neither Uv = —Pq 
nor Uv = —9Pq. Notice that all of these conditions are compatible with our physically 
interesting situation. For —UvPq > is compatible with our choice of signs supporting a 
pair of BPS configurations, and in order to talk about a DODA — D2D6 split none of the 
four charges can be zero. Moreover, since the first term of Eq. (!86|l is positive, the second 
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one for BPS configurations is positive as well {UP > 0), in order to satisfy this condition 
{Uv + Pq){Uv + 9Pq) has to be positive. Hence we see that all of the physically relevant 
conditions are encoded into the structure of the hyperdeterminant 1)4. In particular D4 > 
and —Pvli < of Eq. flM]) gives a necessary condition for the consistency condition to hold. 
Unfortunately the second of our conditions is featuring —Pv which is not coming from any 
of our four-qubit invariants. 

In order to eliminate this shortcoming let us now take another look at the form of our 
polynomial of Eq. (125|) . By employing the substitution 

t = h + y (87) 

this polynomial can be transformed to the reduced form 

y'^ + ay^ + by + c (88) 

where 

a = 6(/2-/i'), 6 = 12/1/2 - 4/3 - S/j', c = -3{I^ - + S (89) 

and S is our well-known quantity defined by Eq. (l29l) or alternatively by Eq. fl68|) . Now the 
conditions 

A4 > 0, a < 0, c < J (90) 

imply^ that our original polynomial Eq.f E^ featuring the fundamental 4-qubit invariants 
is having only real roots. Note that here A4 is the discriminant of Eq. fl25|) . and we also 
recall that A4 = 256/^4. A calculation for the degenerate cases {I4, = 0) shows that these 
conditions are 

> 0, I'l > M, hh > M{ll - M) (91) 
yielding the st^ model and the further specialization M = results in 

L'4 > /1/3 > (92) 

corresponding to the model. Now by virtue of Eq. fl73l) it is easy to see that a further 
specification to the case of our BPS DODA — D2D6 split renders our condition —Pvli < 
equivalent to the one /1/3 > 0. Hence we obtained the nice result that for DODA — D2D6 
splits the conditions encapsulated in the positivity of three-qubit (i.e. U-duality) invariants 
of Eqs.( !7T]) -( l72|) and the positivity of the four-qubit ones of Eq. (p2|) provide a necessary 
condition for the consistency condition for such two-center composites to hold. 
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Let us also verify explicitly that the aforementioned criteria indeed provide real roots of 
our polynomial fl25|) featuring the algebraically independent four-qubit invariants. For the 

model one of the roots is zero due to the vanishing of the invariant J4. For the remaining 
three roots we have to look at the solutions of Eq. fl71|) . After the substitution 

. = t-2(^) (93) 



/^ = -3(^V, ^ = 2^^) -4/3 (94) 



and the definitions 
this cubic equation and its discriminant takes the form 



s-^ + ^s + u = 0, A= (^) + (-) (95) 



where for the explicit form of A see Eq. fjTSj) . According to Eq. f[76l) A < 0, this yields for 
Cardano's formula the case of "casus irreducibilis"— with explicit solutions s^+i, j = 0, 1, 2. 
Transforming back these solutions by using Eq.( p3|) to the variables tj+i we obtain the final 
solutions 

^ ^ (^) (j ^ ' = °' ^' ^ ^^^^ 

where 

a quantity clearly positive by virtue of Eqs. (!73|) i.e. /1/3 > for all sign combinations 
compatible with Eq. flS^ . Since according to Eqs.f lT^ and fITH]) the sign of is just the 
opposite of the sign of A we see that in the special case of the BPS DODA — D2D6 split in 
the model the real roots we have obtained are in accord with our conditions of Eq. (p2|) 
used in a more general context. 



C. Real roots and canonical forms 



Our explicit formulas for the real roots of the fundamental polynomial of Eq. (125|) enables 
an explicit construction of the canonical forms of the four qubit states associated to the 
charge configurations describing two-center solutions. Let us give just a few examples for 
the BPS DODA — D2D6 split. First we write the roots of the fundamental polynomial of 
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Eq.(l25]) for the BPS DODA - D2D6 splits of the model in the following form 

= 4(e - /) sin^ + , J =0,1,2 U = (98) 

e(e + 3/)^ + /(/ + 3e)^ Uv 
e(e + 3/)^ - /(/ + 3e)^ 3 
Note that for BPS solutions we have ef < 0. Notice also that for our special case 2Ji = 

(Qi, Q2) = 3(e — /) cannot be zero for BPS splits so this charge configuration is mutually 

nonlocal. 

As our first example let us consider the nontrivial cases / 7^ 0, and e 7^ when D4 = 0. 
The first case is characterized by the constraint e + 3/ = 0. In this case Is 7^ and from 
Eq.( p9|) we get $ = vr hence 

«,=i3 = ^, i,=4(?|l), i, = 0, A = i|^. (100) 

Using now Eq.( l35|) up to permutations one can associate a canonical form to this configura- 
tion as shown in Eq. (!32l) . This situation arises for example when U = —v = 3 and P = q = 1. 
Notice that this is a highly degenerate case since now D^^Q) = 0, and W = 3(1 — r)(l + r^), 
hence Z = for r = ±1 which is on the boundary of the lower half plane. One can also see 
that for e + 3/ = the superpotential W can be written in a factorized form. 

Our second example is associated with the case when / + 3e = producing the other 
nontrivial zero for D4. In this case Is = 0, and $ = hence we have 

ti = 0, t2 = h = 3('^] , U = 0, h = 2Uv. (101) 



^ 3 , 

From Eq.( l32i) one can see that the canonical state is proportional to the one |0101) + |1010). 

Our last example is a one with $ = |. In this case e(e + 3/)^ + /(/ + 3e)^ = which is 
of the form 

,3,5A^,|,^J_.0. A^^. (102) 

It is easy to check that A = — 1/3 i.e. e+/ = is a solution, hence (A + l/3)(A^ + 14A/3 + l/9) 
is a factorized form of our polynomial, yielding the solutions 

Ai = -^, A2,3 = ^(-7±4v^). (103) 
Let us consider the rational solution A = — |. In this case we get 

ti = ^, t2,3 = |(lTv^)', ^4 = 0, h = Uv. (104) 
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Now the real roots are all different hence ^ 0. However, now Uv = —3Pq hence 
{Uv + Pq){Uv + 9Pg) < 0. A consequence of this is that the consistency condition of 
Eg. (1861) cannot be satisfied so no charge configuration of this kind supports a two-center 
solution. The canonical form can again be read off from Eqs. fl5^ - fl5^ . Notice also that for 
all three cases ti + ^2 + ^3 + ^4 = 4/i as it has to be, moreover the corresponding values for 
A can be written in the form = — l/S'^ for = 0, 1, 2. 

It is interesting to realize that all these special charge configurations giving rise to special 
four-qubit canonical forms are outside the domain of legitimite two-center solutions. Later 
when we connect the special values of e and / to properties of the j function we have some- 
thing more to say about this phenomenon. In order to get charge configurations supporting 
BPS two center solutions we have to chose the asymptotic moduli from the eligible region 
boundend by the usual wall of marginal stability^. In our case this wall is given by the locus 

D. Splitting of invariants 

Notice that the quartic invariant for the model has the explicit form 

- Di{Qi + Q2) = -iUvf + 3{Pqf - GUvPq + AUP^ - 4vq\ (106) 
Now 1^2 = 4:11 P^ and /_i = —Avq^ hence according to Eq.f lM|l we can write 

- /^3(Ql + Q2) = n4(A.on.2i3, 1) = [4/+1 + 6/0 + 4/_i] + 1+2 + 1-2. (107) 

From this expression we see that the relationship between the quantities —D^i^Qi + Q2), 
—P>3{Qi,2) is governed by the combination 4J+i + 6/0 + 4/_i of 5*^(2)0 covariants. Such 
relationships are needed for studying situations when the two-center solution is BPS, but 
the corresponding single center one is not. 

For DODA — D2D6 splits the quantity above can also be written in two equivalent forms 
featuring quantities related to factors of SL{2)q invariants 

- D;{Q^ + Q2) = -SUvPq - 2h{Uv + Pq) - D-^iQ^) - ^3(22), (108) 



or 



- Ds{Q, + Q2) = -2A{Pqf - 2h{Uv + 9Pg) - ^3(21) - D^{Q2). (109) 
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Hence if both of the constituents are BPS i.e. —D^i^Qi) > and —-03(22) > 0, the 
conditions for the single centered system to be BPS too i.e. —D^^Qi + Q2) > are again 
governed by a four qubit invariant Ji or by the factors of another four qubit one namely 1)4. 
It is amusing to realize that all of the factors of are appearing in these two equations. 

As an example one can see that if Ji < and D > and moreover Uv + Pq > (which 
is a factor of D4), then two BPS configurations again yield a BPS one. In the example of 
Bates and Dene^ U = -v = A and P = q = 1, -DsiQi) > and -^3(22) > 0, hence 
1)4 > and Ji < but Uv + Pq < hence --D3(Qi + Q2) = -125 < 0. In this case a single 
center BPS solution does not exist, although the corresponding two-center one does. 



E. The j invariant 

It should be ovious by now that the basic mathematical object giving rise to our observa- 
tions is the fundamental polynomial of Eq. ( l25l) featuring all of our algebraically independent 
four-qubit invariants. For no matter what kind of two center charge configuration we have 
we can construct this quartic polynomial. As a next step we can calculate its resolvent qubic 
of Eq.f l67p . After the substitution u = 2x this polynomial takes the form 2{Ax^ — Sx — T). 
As a next step to this resolvent cubic we associate an elliptic curve of Weierstrass canonical 
form as 

1/2 = 4x^ -Sx-T. (110) 
The j invariant of this curve is defined as 



In particular for the special case of the model we get 



. _ (2/0^ ((2/,)3 _ 24/3)^ 

^ /| ((2/1)3-27/3) • ^''^^ 



After recalling the definitions of Eq. fl99|) a further specification to our case of the BPS 
DODA - D2D6 split one obtains 

17285^ = (2/i)^((2/i)3 - 24/3)^ = 27(e - /)^[27e(3/ + e)^ - 3/(3e + /)Y, (113) 

256/^4 = S^- 27T^ = -27[e(3/ + e)^] [/(3e + /) (114) 
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Now 

27e(3/ + ef - 3/(3e + ff = (3e - 6/)(3e - 6/)(3e - ^s/) (115) 

where 

^■^^ = -5 + 4v^a;^ = e2"^■/^ j = 0, 1, 2. (116) 

Interestingly the number i^i can also be written in the form 

6^-5 + 4^^3(1^) (117) 

moreover, one can check that ^1 + 6 + ^3 = -15, ^16 + 6^3 + 6^3 = 75 and ^i^^a = 3. 
Finally one obtains for the j invariant the expression 

^ = rr^^f— W^^r^n^[(3e - ei/)(3e - 6/)(3e - ^3/)]^ (118) 
[e(3/ + e)^][/(3e + /)^]^ 

Clearly the denominator of j is zero precisely when is vanishing. According to the 
results of the previous subsections for BPS (e/ < 0) composites when this happens the 
consistency condition cannot be satisfied. The numerator on the other hand can vanish 
either when 2/i = (Qi, Q2) = i.e. the charge configuration is local, or when the remaining 
real factor vanishes i.e. when Uv = C,iPq- Since according to Eq.f lll7p .^1 is positive, for 
BPS composites the latter two conditions cannot be satisfied. Hence the j function is not 
having any pathological behavior for the physically legitimate cases of BPS DODA — D2D6 
composites. 

Notice however, that the j function can be made to vanish for non-BPS DODA — D2D6 
composites. In this case we have ef > 0, and one can try to combine the elementary non- 
BPS solutions obtained for the DODA and D2D6 systems^. For such non-BPS composites 
with local charge configurations a generalization of the consistency condition has recently 
been given^-. 

V. ELLIPTIC CURVES 

A. Elliptic curves and four-qubit states 

Let us try to justify that our association of an elliptic curve to the resolvent cubic Eq. (!67|) 
of our fundamental polynomial of Eq.f l25|) is a natural one. For this it is useful to regard our 
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three-qubit charge states \ Qi) and IQ2) with amphtudes given by Eqs. ( 136|) . ( HT]) - (H2 | ) as the 
ones embedded into more general unnormahzed four-qubit ones with complex amphtudes. 
(We win have something more to say about this embedding later.) Hence we can use the 
algebraically closed field of complex numbers to put our considerations in a more general 
setting. 

Let us now notice that according to Luque and Thibon^^ we can write our key quantity 
256D4 = — 27T^ in yet another form by rewriting the polynomial invariants S and T as 

125 = - 2V, 216T = U^ - WV + 21^0"^ (119) 

where 

U = + A{M - L), V = 12{HD-2LM). (120) 

The algebraically independent invariants^ H, L, M, D showing up in these expressions are 
related to our set Ji, J2, /s, I4 aa^ 

h = \H, I2 = ^{H'' + 2L + AM), h = D + ^HL, h = L. (121) 

It is known that a nonsingular projective curve of genus 1 is isomorphic to a plane cubic 
curve of the Tate fornt^ 

X0X2 + aiXoXiX2 + asXQX2 — x^ — a2XoX-^ — a^x^Xi — UqXq = 0. (122) 

Here the coefficients 01,02,03,04,06 can be taken from an algebraically closed field F of 
arbitrary characteristic and (xq, xi, X2) are homogeneous coordinates of the projective plane 
FP^. (Of course for qubits we are merely interested in the complex case i.e. our concern 
will be curves in CP^.) Using inhomogeneous coordinates {x,y) = (xi/xq, a;2/a;o) this can 
be rewritten as 

y"^ + aixy + a^y = x^ + 02X^ + 04X + 05. (123) 
Now the j invariant and discriminant A of a plane cubic curve of this form is defined^ as 

J = I = 1728^^ (124) 

where 

C4 = 62 - 2464, C6 = -&2 + 3662^4 - 21666 (125) 

&2 = ai + 4o2, 64 = 0103 + 204, 60 = c^s + 4o6. (126) 
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A necessary and sufficient condition for our curve to be nonsingular is A 7^ 0. 

Now looking at Eqs. flll9p - fll20p one can easily check that four-qubit states provide a 
particularly nice parametrization for a special class of plane cubic curves with qq = 0. 
Indeed setting 

b2 = U = + 4{M - L), b4 = = HD -2LM, be = (127) 

yields 

C4 = 12^, C6 = -216T. (128) 

Moreover, one can also see that 

ai = H = 2Ii, a2 = M-L, a3 = D = I3- hh, = -LM, = (129) 

hence the family of such cubic curves for the stu model {H ^ 0, D 0, M ^ 0, L ^ 0) is 

2/^ + Hxy + Dy = x^ + {M - L)x^ - LMx (130) 

for the st^ model (i? 7^ 0, L> ^ 0, M ^ 0, L = 0) 

y'^ + Hxy + Dy = x^ + Mx"^ (131) 

and finally for the model (// 7^ 0, D 7^ 0, L = 0, M = 0) it is 

1/^ + Hxy + Dy = (132) 

For all of these curves the hyperdeterminant is given by the usual formula 256D4 = 
— 27T^ with the j invariant having the form of Eq. flllip . As an example let us take a 
look at the hyperdeterminant D4 for the st"^ model which is given by the expression 

256^4 = /|(/3[(2/i)^ -27/3] +36M[/2 -2/1/3]) , 2M = 3/2 - 2/^. (133) 

One can see that in accordance with our association of cubic curves to the stu model and 
its truncations for M = we get back to Eq. (176|l i.e. the expression for the model. 

Recall also that if the characteristic of the field is not 2 or 3 via further projective 
transformations (i.e. completing the square on the left hand side and then the cube on 
the right hand side of Eq. (11301) ) one obtains the Weierstrass canonical form of Eq.f lllOp 
which is simply related to the resolvent cubic of our polynomial S4 of Eq. (l25p we have 
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started our considerations with. Notice however, that although due to the fact that now 
F = C the Weierstrass form in our case can always be reached, this form is not showing the 
basic differences between the structures of the cubic curves associated to the stu, st^ and 
models, hence we prefer the Tate form. 

In the following as usual by an elliptic curve we will mean a nonsingular plane cubic 
curve as given by Eq. fll22p and the "point at infinity" (xo,xi,X2) = (0,0,1). Since for 
these curves we have only (0, 0, 1) as a point at infinity by an abuse of notation one can 
call Eq. fll23p as the defining equation for the corresponding elliptic curve. An elliptic curve 
is defined over a subfield F' of F if all the coefficients 01,02,03,04,06 are taken from F' . 
In the case of the discrete version of the continuous [/-duality group, as nonperturbative 
string symmetries, the coefficients of this curve are the algebraically independent four-qubit 
invariants if, L, M, D (related to the ones Ji, J2, /a, 1^ via Eq. fll2ip ) taken from the ring of 
integers of the subfield Q. In this case the arithmetic aspects of the theory of elliptic curves 
within the context of two-center black hole solutions becomes important^. 

The upshot of these considerations is that to a particular charge configuration of two- 
center black holes of the stu model we can indeed associate quite naturally an elliptic curve 
of special kind (i.e. Og = 0). This elliptic curve is of the form of Eq.f ll30p . 

Notice that the invariants H, L, M play an important physical role. Since H = {Qi, Q2) 
the vanishing of this invariant gives rise to local charge configurations, in this case the xy 
term of Eq. fll30p is missing. The vanishing conditions on L and M give rise to truncations 
to the st"^ and models. Notice also that since for L = we have D = hence in the 
st'^ and models sending D to zero corresponds to taking the limit when our elliptic curve 
degenerates (see Eqs.([76D, (fT2T|) and ffTSSj) ). 

In establishing the correspondence between two center charge configurations and elliptic 
curves we tacitly assumed that our unnormalized charge states with real amplitudes (trans- 
forming according to the corresponding representation of the continuous [/-duality group 
valid in the supergravity approximation) are regarded as ones embedded in unnormalized 
four-qubit states with complex amplitudes. However, until this point we did not specify 
what kind of complex states we are having in mind. 

In a recent series of papers it was shown^^i^ that for single center extremal stationary 
spherically symmetric BPS and non-BPS black hole solutions apart from charge states i.e. 
unnormalized three-qubit states with real amplitudes, we can also introduce unnormalized 
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complex thiee-qahit states which are depending on the charges and also on the moduh fields. 
(For a more general class of such states even the warp factor can be included^.) It has been 
shown that the complex amplitudes '?/'ooo ; V^ooi > • • • > V^iii of such states {ip) are related to the 
central charge and its covariant derivatives with respect to the moduli. For example the 
amplitudes featuring the central charge have the form^^ 

^000 = V2Z, ^111 = -V2Z. (134) 

Moreover, recently it has also been demonstrated that such complex three-qubit states are 
just special types of four-qubit ones^^. It is natural to expect that for two center solutions 
complex four-qubit states of similar type should show up. For example in the consistency 
condition of Eq.( l77|) the term Im(2^22i) can be reinterpreted as a truncation of a pairing 
between the two complex amplitudes \I'oooo = V^^i and = —\/2Z2 of such four- 

qubit states. In this respect recall the formalism developed by Denef^ in the type IIB 
duality frame where apart from the usual antisymmetric, topological moduli independent 
intersection product (our four-qubit invariant Ji in the stu model) the importance of the 
symmetric, positive definite, moduli dependent Hodge product is emphasized (in the stu 
model giving rise to the interpretation of the black hole potential as the norm of a three- 
qubit stated) . Such structures are the natural ones appearing in the multicenter black hole 
context. In the stu case it is easy to establish a correspondence between this formalism 
and the one where complex three and four-qubit states appear—. We conjecture that such 
moduli dependent complex states and their associated invariants featuring elliptic curves 
could be the relevant mathematical objects hiding behind the considerations of this paper 
based merely on charge states. Note that moduli dependent curves of that kind would also 
display explicit dependence on extra complex parameters giving rise to elliptic fibrations. 
We will have something more to say about this possibility in Section V. 

B. Degeneracies 

The discriminant A of our elliptic curve is related to the hyperdeterminant as 2560^ = 
A. For nonvanishing the curves are nonsingular. Since our curves are of genus 1 they are 
topologically tori. Vanishing of 1)4 results in different degeneracies of these tori (for example 
one of their homologically nontrivial cycles can contract to a point). One can illustrate this 
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in the stu and st"^ models when both of the two BPS constituents are having vanishing 
Bekenstein-Hawking entropy i.e I±2 = for both centers. 

Indeed, let us use the parametrization as introduced by Sen^ for such two-center compos- 
ites Q = Qi + Q2 



adQ -abP\ / abP - bcQ \ ah 
cdQ - bcP j \adP -cdQj ' \cd 
Here P and Q are four- vectors related to the ones of Eq. (j38l) -(l42ll as follows 



eSL{2,Z). (135) 



A = a{dQ -bP), B = c{dQ-bP), C = b{aP -cQ), D = d{aP - cQ). (136) 

Since the vectors A and B and C and D are now proportional the bivectors of X and Y of 
Eq. fH3l) are zero hence I±2 = -^±1 = though Iq 7^ 0. Now Eq. fl68|) shows that 

S = Ml T = II (137) 

hence — 27T^ = i.e. D4 = 0. The nonseparable bivector is Z = A P hence 

h = lz' = ^{QAP)-{QAP) = -^D.iQ). (138) 

Alternatively one can directly check that /a = J4 = and 6/2 = —6/0 = D^IQ) in accordance 
with Eq.f l56|) . hence S = 3/| and T = /| yielding again D4 = 0. 

In order to obtain the Tate form of our elliptic curve we have to calculate the quantities 
H, L, M, D. One can show that 

H = P-Q, L = 0, M = ~P^Q\ D = (139) 



where we introduced the notation 

'q\ _ ( d -l> 
P) - [-C a 

Now the elliptic curve of Eq. (11301) is of the form 




(140) 



y^ + P- Qxy = x^- ^P^Q^x'^. (141) 

Completing the square on the left hand side (and by an abuse of notation using y for the 
new variable again) yields 

.221 ^ 



y' = x'[x + -D,iQ)]. (142) 
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Notice that the quantities P and Q hence the Tate form of the cubic curve is featuring 
the parameters a, 6, c, d characterizing the two-center spht, however the latter form is not. 
This is due to the fact that D^{Q) = {P ■ Qf - P^Q^ = {P ■ Qf - P^Q\ This leads 
us to the important observation that at least for this singular example two-center charge 
configurations labelled by elements of SL{2, Z) correspond to the same singular curve with 
the (I142p canonical form . Hence apart from the four copies of 5'L(2)s related to four-qubit 
systems there seems to be a fifth SL{2) parametrizing a hidden torus. We will return to this 
important point later. 

Eq. fll42p represents a cubic curve with a node at the point (0,0). The tangent lines at 
(0, 0) are obtained from the equation y"^ — [D-^/^^x^ = 0. Such lines are of the form 

y-l^/D^] (y+ly^] =0. (143) 



2 ^ 7 V 2 

It is important to recall that we assumed that all of our curves are over the complex 
numbers, i.e. F = C. Hence for a cubic in the canonical form y"^ = + ax + (3 we are still 
free to use transformations of the form {x,y) h^- (/i^x,/i^?/) and (a, /3) (-)■ (fi'^a, fi^f3), where 
/i G C^. Using the transformation {x,y) (— x, —iy), i.e. fi = i transforms Eq. (11421) to the 
form 

y' = x'(x-]D,{Q)]. (144) 



4 

As we see such transformations are similar to the ones changing a BPS charge config- 
uration {D3 < 0) to a non-BPS one {D^ > 0). Interestingly as shown by Eq. (ll43p (and 
a similar one obtained from Eq.f ll44p ) the pair of lines at the node are parametrized by 
a/— D3 for the BPS-case or a/I^ for the non-BPS one, i.e. quantities proportional to the 
Bekenstein-Hawking entropy. 

Notice also that Eq.( ll44p can also be parametrized as 

y^ = 3t{x'-tf + {x'-t)^ = {x')^-3t\x')-2t\ ^=-^ (145) 
Using Eqs. fll37p - (ll38p we can write this as 

y' = x'' - \sx' - It. (146) 

Now a further transformation {u,v) = {x/2,x/8) i.e. /i^ = 1/2 transforms this to the form 
v"^ = — Su — 2T of Eq.f l67p . with the left hand side being just the resolvent cubic of 



32 



S4 = 0. The procedure discussed here illustrates how the general Tate form of our cubic 
curve of Eq. (11411) can be transformed to a resolvent cubic form. 

It is important to realize that if the domain of definition of our curve is the subfield Q 
from the general Tate form of Eq. (11301) merely the form 

y^ = x^-\sx+^-T (147) 

can be reached. Further reduction is possible only if /i is a square in Q. In the case of our 
singular curve of Eq. (I142p if is not a square we cannot represent the tangent lines in the 
form of Eq. (ll43p hence the tangent lines in this case are not rational. 

Due to permutation symmetry we have no parametrization of (11351) type for the model 
available. However, as a degenerate example in this case we can consider the example of 
Bates and Denef^ instead by sending one of the charges e.g. P to zero. In this case one of 
the black holes is a small one, and the two center composite is characterized by the charges 
U,v,q. Now one checks that I3 = 0, hence according to Eq.( 176i) we get again D4 = 0. The 
explicit form of the corresponding degenerate elliptic curve is 

2/2 + Hxy = x\ H= (Qi, Q2) = Uv, (148) 

with j = 00 just like in the previous example. 

C. Invariance properties 

Since the coefficients aj of the elliptic curve of Eq. (11301) are invariants the same curve is 
associated to equivalent two-center charge configurations in the stu model. In other words 
our curve is clearly invariant under the action of the group S'L(2,R)o x G4, where SL{2,M.) 
is the horizontal symmetry group of generalized exchange transformations of the centers and 
G4 is the d = 4 continuous ?7-duality group of the stu model which is SL{2,W)^^. 

However, as far as physics is concerned in a special case of the model we have also 
found some connection between the structure of SL{2,M.)q x 6*4 (four-qubit) invariants and 
the structure of the consistency condition of Eq. (l77l) . This connection at first sight is not 
surprising since the consistency condition contains a four-qubit invariant 2/i = {Qi, Q2), 
however more importantly it is also featuring the quantity Im(22^i) which is not invariant 
with respect to the full group of 5'L(2,M)o x G4 transformations. 
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In order to show this let us observe however, that this quantity is invariant under trans- 
formations belonging to the horizontal subgroup SL{2, M.)q. For let us recall our definition 
of Eq. and then write the central charges as 

Zi = e^/2(Qo + Qm + Q2r2 + QsTs - P^t^t^ - P^nr-, - PViTs + P'nr^rs) (149) 

Z2 = e^/^(go + + g2T"2 + - ^^^"2X3 - pViTg - p^TiT2 + P^TiT2T'i), (150) 

then for a,b,c,d E M we have Z[ = aZi + bZ2, and Z'2 = cZi + dZ2, hence lm(Z'2Z[) = 
lm{Z2Zi) due to ad — be = 1. However, this quantity is invariant merely under a special 
subgroup^^ of G4. Indeed, we have 

Z{Kn, Kt2, Kts; KQi) = Z{n, r2, r-,- Q,) (151) 

where K is the subgroup of transformations of the form 

Ta^Ta + K K^R, a = 1, 2, 3. (152) 

Notice that the group K with ka E 'L is precisely the stabilizer of the cusps in the three 
copies of the fundamental domain of the modular group, i.e. the stabilizer of Tq = —ioo for 
a = 1,2,3. 

Since the walls of marginal stability are determined by the central charges and for the four 
qubit invariant 2/i we have {KQi, KQ2) = (Qi, Q2) clearly wall-crossing does not obstruct 
the K subgroup. Therefore we see that unlike our elliptic curves the consistency condition 
is left invariant (in the above sense) merely with respect to the subgroup SL{2,M.)o x K. 
Hence the correspondence between the physics of two center solutions and our special class 
of elliptic curves should be refined (see in this respect the comments at the end of Section 
V.A). 

Now we turn to another important issue we have not discussed yet. Based on our ex- 
perience with the degenerate case studied in the previous subsection we expect that for 
the nondegenerate case our mapping of two center charge configurations characterized by 
four-qubit invariants to elliptic curves should be many to one. Moreover, since in the Tate 
form we should have = our mapping in the stu case cannot be onto either. We also 
know from the theory of elliptic functions that if we work in the algebraic closure of our 
field F i.e. F^^^ elliptic curves are isomorphic if and only if their j invariants are the 
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same. Hence in the case of M working with the complex numbers charge configurations with 
different four-qubit invariants could be mapped to isomorphic elliptic curves with the same 
j invariant. In the singular case according to Eqs. fll39p - (11401) the values of H and M are 
clearly different for different splits, however j = oo in all cases. 

Now over C every elliptic curve E is isomorphic to an elliptic curve E{A) where A is 
a lattice in the complex plane and the mapping between C/A (a torus) and E{A) is an 
isomorphism provided by the usual map defined by the Weierstrass V function. Explicitly 
E{A) is of the form 

xqxI = Axl - g2if)xixl - g3iT)xl (153) 
where 5'2(^) = 60G'4(f) and gsi^f) = lAOG^^f) with G'2fc(^) are the Eisenstein series 

and if X = xi/xq and y = X2/X0 then x = V{z) and y = V'{z). Here the lattice vectors of A 
in C are wi = 1 and uj2 = f with f G EI being the modular parameter of the torus and the 
Weierstrass V function is defined as 

1 1 



[z + m + nfy (m + nfY 



(155) 

I Z -I- 711, -+- ti l \- I III. -+- III \- 

(m,n)GZ2,(m,n)^(0,0) 

Here we have used the notation f for the modular parameter in order not to confuse it with 
the moduli ri, r2, of the stu and with the moduli r of the model. Now the j invariant 
gives rise to the j function 

^■(^) = 1728 3^^^;j! (156) 

Comparing this with the expression of the j invariant as given in terms of the four-qubit 
invariants S and T of Eq.l lllip to a two center charge configuration we can associate a torus 
with modular parameter f. Now j : HI — t- C is an automorphic function i.e. 



cf + d 



a h 



Hzr—A=3{r), gSL(2,Z). (157) 



c d 



It can be regarded?^ as a map providing a holomorphic universal covering of the orbifold 
'H/PSL{2,Z) (i.e. the fundamental domain for PSL{2,'Z)). Its Fourier series at f = ioo is 
provided by the variable u = exp 27iif as 

j{f) - 744 = - + E c„m", c„ G Z+ U {0}. (158) 



n=l 
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From this we see that the two-center charge configurations with discriminant D4 = of 
the previous subsection with j = 00 should correspond to the modular parameter f = ioo 
i.e. the cusp of the Riemann surface 'H/PSL{2,'Z). The different possible two-center splits 
labelled by elements of SL{2, Z) not changing the value of the j invariant seem to be related 
to the modular property of the j-function as shown by Eq. fll57p . It is also intriguing to 
recall that the tangent lines at the node of the degenerate elliptic curve are parametrized by 
the Bekenstein-Hawking entropy (see Eq. fll43p ). 

From these investigations it is natural to conjecture that similar invariance properties of 
other two-center splits with D4 ^ should hold. It would be nice to uncover the physical role 
of this hidden torus and the extra SL{2) symmetry associated with it. Some speculations 
on how the extra modular parameter f should be implemented into a four-qubit picture will 
be given in Section VI. 

D. A conjecture 

Let us consider the elliptic curve of Eq.f ll30p associated to two-center charge configura- 
tions of the stu model. As we have noticed this curve is of the Tate form satisfying the 
special constraint ag = 0. The coefficients aj are four-qubit polynomials with definite degree 
of homogeneity 2j. Indeed, oi = if is a quadratic, 03 = D is a sextic, 02 = M — L is a 
quartic, and 04 = —LM is an octic polynomial in the amplitudes of the four-qubit charge 
state. (The variables x and y should be assigned the degrees four and six respectively.) Now 
according to this observation the constraint ag = should be arising from the vanishing 
condition of a polynomial of degree 12. It is easy to identify this polynomial constraint. It 
is just the constraint^ V12 = of Eq.( l57|) . or alternatively of Eq.( l58|) . hence in a redundant 
notation for the stu model we can write 

+ Hxy + Dy = x^ + {M- L)x^ - LMx + V12. (159) 

Since the stu model can be regarded as a consistent truncation of = 8, d = 4 maxi- 
mal supergravity with continuous^^ [/-duality group G4 = -£^7(7) one might conjecture that 
it should be possible to substantially generalize our considerations concerning two-center 
charge configurations by studying elliptic curves of the form 

y^ + P2yx + P^y = + P^x^ + P^x + P12 (160) 
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where P2j are polynomial invariants of the group SL{2)q x G4 where SL{2)q refers to the 
horizontal symmetry groupi acting as a generalized exchange symmetry group on the two 
centers. Now we have two sets of 56 component charge vectors where each of these vectors 
is transforming according to the symplectic irreducible representation of 6*4 which is now 
the fundamental of -£'7(7). According to an analysis of two centered magical charge orbits'^ 
working in the complexification of G4 one discovers that the two centered charge orbits 
correspond to different real forms of the quotient of the complex groups Ej/ S0{8). In par- 
ticular we have two |-BPS two-centered charge orbits with one of them being Ej^y^/ SO{A, 4). 
Recall now Eq.([l]) displaying in its structure the group theoretical reason for our occurrence 
of four-qubit states. According to this equation one should be able to obtain our four-qubit 
invariants as special cases of SL{2)q x G4 ones. Hence we conjecture that an elliptic curve 
of the f ll60p form should display this reduction procedure via implementing the constraint 
V12 = via a suitable reduction of some unknown SL{2)q x polynomial invariant P12. 

In fact many pieces of such invariants are already at our disposal. For example the 
N = 8 analogues Jo, J±i, J±2 of the four-qubit covariants Jq, /±i, I±2 of Eq.(H9l) arising from 
the polarization of Cayley's hyperdeterminant can now be obtained from the polarization of 
Cartan's quartic invariantr^. Using these explicit expressions one can construct the quanti- 
ties answering the similar ones of Eq. (168|) 

S = 3Jo^-4J+iJ_i + J+2J-2, r = J^ + JlrJ-2 + J^-iJ+2-J+2JoJ-2-2J+lJoJ-l. (161) 

Indeed precisely these polynomials of order 8 and 12 has been suggested as obvious candi- 
dates for members of a complete basis for SL{2)q x G^ invariants. Now one can define an 
elliptic curve 

y' = x^- ^Sx + -J (162) 

as the one corresponding to Eq. fll47p . As was commented there this form can always be 
obtained from the Tate form if the characteristic of the field is neither 2 nor 3. Hence we 
conclude that our Tate form Eq. fll60p featuring the unknown invariants P2j should reduce 
to Eq. fll62p after completing the square on the left and the cube on the right hand side. 
Explicitly this process amounts to the transition from using the coefficients aj = P2j to using 
the ones C2 = - jS and C4 = of Eqs.(^2B-UM- 

As a solid piece of evidence it is also obvious that P2 = {Qi, Q2) with (, ) being the usual 
symplectic product of charge vectors which is a singlet with respect to SL{2)q x G4. Hence 
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for mutually local charge configurations our cubic curve of Eg. (11601) falls short of the term 
proportional to xy just like in the stu truncation. Moreover, as shown by Eqs. (ll2p -f lT^ in 
the four-qubit case our quadratic and sextic invariants Ji and are duals of each other 
. Hence we expect that the invariant Pg should be related to a quadratic combination of 
the Freudenthal duals of the corresponding charge vectors which are cubic in terms of the 
original charges^. Moreover, these dual quantities should be antisymmetric with respect to 
the exchange of the centers. Luckily a quantity Jg? satisfying these criteria is also available 
(see Eq. (3.23) of Andrianopoli et.al.— ). 

Unfortunately since according to Eq.f ll2ip = 1^ — /1J4 in the stu case this quantity is 
probably not directly related to our invariant Pq. Hence one is still left with the problem of 
finding the invariants P4, Pg and P12. 

In order to gain some insight into the physical meaning of the invariant P4 let us look 
again at the stu (four-qubit) case. In this case we know that the corresponding invariant 
is M — L. As a first step we would like to somehow relate this quantity to some invariant 
known within the context of = 8, ci = 4 supergravity. 

For the stu truncation let us suppose that we have chosen a particular two-center charge 
configuration parametrized by the four real numbers a, 6, c, d of the canonical form of Eq.f p2|) . 
Now the explicit form of the invariants H, M, L, D is 



H=t^a' + b' + c' + d') 



M 



c — d 



c + d 



(163) 



D = -{ad — bc){cd — ab){ac — bd). 



abed. 



(164) 



Consider now the quartic -£^7(7) invariant^ expressed in terms of the matrix of central 
charge 2^ of A^ = 8 supergravity in d = A. This matrix is an 8 x 8 complex antisymmetric 
one which can be brought to the canonical form after using a suitable U (8) transformation 
Z H> U'^ZU. This canonical form is 

f zo 0\ 
^1 
^2 
yo 23/ 



U^ZU 




(165) 
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Now Cartan's quartic invariant is 

JiiZ) = Tt{ZZZZ) - ^{TT{ZZ)y + 4(Pf (Z) + Pf(^)) (166) 

where the Pfaffian is 

Pf (Z) = J^^^^'^^^^^^Z^B^CD^EF^GH, A, 5, ■ ■ ■ = 1, 2, . . . 8 (167) 

Subscripts A,B... label an of 8 SU{8). Using the canonical form one gets the well-known 
expression 

Ja = (kor + kiT + + ksH^ - 4:{J2 + 8Re{zoZiZ2Z3) (168) 

It is also known that by an SU (8) transformation it is possible to remove three phases, so 
for instance 2:1,2:2,^3 can be chosen to be real (or else having the same phase). Hence the 
canonical form can be characterized by four real numbers r„ = \zn\ and a phase 0. The form 
of J4 reflecting these considerations is 

J4 = [('^o + ri)^ - + r^f] [(ro - ri)^ - (ri - rg)^] + 8rorir2r3(cos 0-1) 
= -16M + 8L(cos - 1) = 8(A^ - M) + 8L cos 

where we have used the identity L + M + N = Oof Eq. fl22l) and we made the identifications 
ro = a,ri = b,r2 = c and = d. Hence though we did not manage to get directly to 
the invariant M — L however for = | we get N — M instead which is related to M — L 
via a permutation of the qubits labelled by ziZ2^3 carrying the G4 labels. Notice also that 
the permutation group S3 responsible for this triality symmetry at the stu model level is 
related to triality of the complex group 5*0(8, C) whose real form 5*0(4,4) is featuring the 
tripartite entanglement of seven qubits interpretation of Cartan's quartic invariant^^i^ of 
N = 8 d = 4 supergravity and the corresponding stu truncation. Notice also that according 
to Eq.( |56l) M — is just the invariant x of Ferrara et.al.-. These considerations show that 
the unknown invariant P4 should truncate to x = ^ ~ after some suitable permutations. 
It is also interesting to notice that for Zn real i.e. = we have 

J4 = -16M = 16forir2r3. (169) 
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Here 



m=0 







/l 1 1 l\ 


/ro\ 


fi 


1 


1-11-1 


ri 




~ 2 


11-1-1 


^2 






V -1-1 1 / 





(170) 



where if is the Hadamard matrix used for implementing discrete Fourier transformation 
in quantum information theory. Notice that with the definition Qn = 2f„ we have J4 = 
Q0Q1Q2Q3 the standard expression giving rise to the Bekenstein-Hawking entropy S = 



'^V^QoQiQ2Q3 of the D2 — D2 — D2 — D6 brane configuration in the type HA duahty 
framed. 

There still remained the task to understand the meaning of Pg and Pu- Since in the 
N = 2 , d = 4 context the vanishing of the corresponding invariants identifies the stu, st^ 
and truncations in terms of the structure of the fll60p elliptic curve it would be important 
to know their explicit forms. Finally the discriminant — 277"^ as a polynomial invariant 
of order 24 should play the role of some sort of generalization of the hyperdeterminant of 
type 2x2x2x2. It would be interesting to clarify what kind of role this discriminant and 
the associated j function plays in the physics of two-centered black hole solutions. 



VI. A TRIALITY SYMMETRIC CURVE 

The aim of this speculative section is to draw the readers attention to some interesting 
structural similarities showing up in a variety of physical contexts where our four-qubit 
invariants parametrizing elliptic curves might play a crucial role. 

In the previous investigations our basic philosophy for the association of elliptic curves 
to two-center charge configurations was based on the resolvent cubic of the fundamental 
polynomial of Eq. fl25p . Is there any other physically appealing way for this association 
which retains the fundamental role of this polynomial and at the same time also featuring 
an extra modular parameter f accounting for a hidden torus? In order to show that the 
answer to this question is yes let us rewrite Eq.f l25p in the form 



S4(A,„,,,,i3, f, -x) = + («/3)(4Ji)x^ + {a(3)\<6h)x' + {al3y{Ah)x + {al3fli. (171) 
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The quantities displaying explicit dependence on f are defined as 



a = 62 - ei, /3 = 63 - ei (172) 

where 

4(x - ei)(x - 62) - eg) = - - 5(3 (173) 
with g2 and given by Eqs. fll53p -f ll54p . Explicitly we havei^i^ 

ei-e2 = ^3^(0,f), e3-e2 = ^t(0,f), - = et{Q,f) (174) 

where 

^,(0,f) = 5^g^("+i/2)2 (175) 

^2(0,f) = 5^(-l)V"' (176) 

neZ 

^3(0,f) = ^g^-^ (177) 



neZ 

with g = e^'^*'^. One also has^"^ 

ei = ^ + 16g + 16g2 + . . . (178) 

62 = -- - 8g^ - 8g - 32gt - 8^^ + . . . (179) 

3 

63 = -- + 8g^ - 8g + 32gi - 8g2 + . . . (180) 

3 

Notice that 61 + 62 + 63 = similarly to the property of four qubit quartic invariants 
L + M + N = Q. In the limit f — )■ ioo we have a/3 — ?■ 1 hence from Eq. fll7ip we get back 
to our usual polynomial of Eq. fl25p . Clearly introducing the factor a/3 amounts to rescaling 
the amplitudes of our four-qubit state of Eq.([2]) from A to the f dependent ones a/o^A. 

Let us now introduce an extra complex parameter u to be specified later and consider 
the following family of elliptic curves 

= x{x — au){x — (3u) — (a — /3)^/ia;^ + -(a — /3)q;/3[(q; + /3)/4 — 3(q: — (5)l2\x 

- {a- P)a^p'^hu- {a- Pfa^p'^h. (181) 
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Now the right hand side is a quadratic polynomial in u. Its discriminant turns out to be 
just (a — /3)^S4(A, f, —x). This gives the desired clue for yet another way of associating an 
elliptic curve to two-center charge configurations. 

Notice that in the limit f — )■ zoo our curve boils down to = x{x — uY a curve similar to 
the one of Eq. f ll42p with j invariant oo. Moreover, via the transformation v = u + {a — f3)^H 
one can also obtain a new form 

y2 = _ (yv){x " /3v) + l3{a - /3)Hx^ - af3{a - f3){aM + f3N)x 

- al3{a - P)Hvx - a^(3\a - P)Lv - a'^P^a - P^D. (182) 

Though not in the Tate form this curve is displaying similar quantities than our previous 
curve of Eq. fll30p . In particular for truncations from stu charge configurations to st"^ and 
ones the structure of the curve is getting simpler step by step. Though this curve now 
shows an explicit dependence on a modular parameter f and also featuring a new complex 
variable v, now its discriminant cannot obviously be related to a hyperdeterminant of type 
2x2x2x2. However, curves like Eq.f ll82p have other appealing properties which we would 
like to discuss. 

Actually the curve of Eq. fllSip is related to a one introduced by Seiberg and Witten 
in their study of = 2 supersymmetric SU{2) gauge theory with four quark flavours^. 
The parameter u in that case was the gauge invariant modulus representing the square 
of the Higgs expectation value and f was the complex coupling constant O/n + Sni/g"^ of 
Montonen- Olive duality. In that context their curve was parametrized not by algebraically 
independent four-qubit invariants but by the squares of the four quark masses mi,m2,ms 
and m^. However, after comparing the relevant expressions (see in particular Eq. 17.58 of 
Ref.— ) it is easy to show that the quark masses squared in that case correspond to the 
squares of the parameters a, b, c and d of our four-qubit canonical form of Eqs. fl32|) - flM|) . 

Apart from the symmetry SL{2,'R)q x G4 where G4 = SL(2,M.)^^ is the continuous U- 
duality group in the supergravity approximation, the stu model also has an important triality 
symmetry^. In our four-qubit description this 5*3 permutation symmetry should manifest 
itself in representing somehow elliptic curves in a way displaying four qubit invariants that 
are also invariant under permutation. In analogy with the Seiberg- Witten curve there is also 
the possibility to present a triality invariant form of that kind. For this purpose we have 
to chose from four algebraically independent four-qubit invariants which are also invariant 
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under the full permutation group 5*4. Such invariants were first constructed by Schlafii^ in 
1852. Let us denote these invariants as^i^ if, F, S and 11. They are of order 2, 6, 8 and 12 
respectively. The new quantities F, E and 11 are defined as 

T = D + E + F, S = L2 + M2 + Ar2, U = {L - M){M - N){N - L). (183) 

Here E and F are defined as^ 

HL = E-D, HM = D-F, HN = F - E. (184) 

Notice in particular that in terms of the quantities H, L, M, D used to describe the Tate 
form of the elliptic curve of Eq. (11301) the permutation invariant version of the sextic invariant 
reads as 

T = 3D + H{L- M). (185) 

Let us now consider the triality symmetric curved 
A 

= yj^W2W3 + -[{M-N){e2-e3)wi + {N-L){e3-e,)w2 + {L-M){ei-e2)w3]-Y^ (186) 
where 

Wi = X - eiU - e^iH, A = (ei - e2)(e2 - e3)(e3 - ei) i = l,2,3. (187) 

Notice that after using Eq. f ll85p and the canonical form for our four qubit states of Eq. ( 132|1 
and fl35|) with (ti, t2, ts, ^4) = (a^, 6^, c^, cP) for the polynomial F of sixth order we obtain the 
expression 

-F = — V Udk - — V Ut^i + — V . (188) 

i>j>k i^j i 

This is precisely the sixth order invariant of Eq. (16.36) of Seiberg and Witten^ provided 
we make the identification ti = i.e. the canonical four-qubit parameters are identified 
with the quark masses of the four fiavours. Similar calculations verify that the remain- 
ing invariants of that paper namely R,Ti,T2 and T3 can be identified with the invariants 
H, |(M — A^), ^{N — L) and ^iL — M) in the canonical parametrization. The important 
property of the curve of Eq. f ll86p is that after making the transformations 

u = u + 2^iH, X — )■ z — + 2^1H (189) 

and performing the (weak coupling^'') limit f — )■ ioo it boils down to the form 

y^ = x'^{x-u). (190) 
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Moreover, in order to make contact with our original curve of Eg. (11811) one just has to 
perform the change of variables in Eg. (11 86 p 

X I— )■ X — eiu — e\H, u = u + ij^iH. (191) 

It is also interesting to realize that our identification of the guark masses sguared with 
the canonical four-gubit parameters automatically incorporates the special cases when some 
of the masses are zero with the singular cases in the four-gubit classification scheme of 
Verstraete^^. Moreover, switching to the two-center stu context there the canonical param- 
eters would rather be identified with the the parameters ^1,^2,-23 and 2:4 of the canonical 
form of the central charge matrix of = 8 ci = 4 supergravity (see also the discussion of the 
previous seubsection) . Recall also in this respect the observation of Seiberg and Witteni^ 
that triality symmetry of S'0(8) is connected to the permutation symmetry 5*3 as the mod 2 
reduction of SL{2,Z). Comparing Egs. (17.34)-(17.35) of Ref.l4 with our Eg. (11701) related 
to the structure of the matrix Ai and a similar expression related to the matrix Af that 
are in turn related to the structure of four-gubit reduced density matrices(see Egs. ([7])- ([8])) 
we see that triality symmetry is intrinsically related not only to the structure of four-gubit 
entanglement but to the structure of the manifold X described by the curve of Eg. (1 186 p . In 
this context it is especially instructive to recall the arguments of Seiberg and Witten on the 
structure of the cohomology of X . In particular it is tempting to reinterpret the triality 
invariant expressions for the cohomology classes of the periods^ as genuine unnormalized 
four-gubit states in canonical form (see Egs. 17.36 and 17.37 of that paper). The redun- 
dancy in arriving at such an expression can be attributed to the action of the Weyl group of 
5*0(8) on the canonical parameters (i.e. the action of the Klein group on the canonical form 
of |A) in accord with a comment of Lugue and Thibon in the four-gubit contexir^. (See also 
the last paragraph of Section II. in this respect.) 

Finally it is amusing to recall yet another context where elliptic curves parametrized by 
four-gubit invariants reveal some intriguing structural similarities with interesting physics. 
First of all it is well-known that in F-theory an extra SL{2, Z) and a hidden torus also 
makes its presence in a spectacular way^^. In this case one considers an elliptic fibration 
M with some basis manifold B and fiber a two dimensional torus with modular parameter 
f. Now F-theory is defined on M as type IIB string theory on B with the axion dilaton 
modulus of type IIB string theory identified with the modular parameter of the two-torus. 
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In the original setting^ an elliptically fibered K?) surface was considered of the form y"^ = 
+ f {u)x+ g{u) where / and g are polynomials of degree eight and degree twelve in u G CP^. 
This curve describes a torus for each point of the Riemann sphere CP^ labelled by the 
complex coordinate u. It is interesting to realize that this elliptically fibered K3 surface is 
of the form of Eq. (11471) where S and T are polynomials of order eight and twelve. However, 
according to Eq. fl32l) when writing S and T in the canonical form we are having four complex 
coordinates instead of the one u. The points where the torus degenerates corresponds to the 
vanishing of the discriminant of the cubic which is a polynomial of order 24. In the four qubit 
parametrization this just corresponds to the vanishing of our hyperdeterminant of order 24 
related to the quantity — 27T^. In the original F-theory context the compactification 
described by this elliptical fibration corresponds to a configuration of 24 seven branes of type 
IIB theory located at the zeros of the discriminant. On the other hand in the four-qubit 
canonical parametrization the explicit expression for the discriminant of Eq. (135!) is similar to 
the usual ones describing coincident seven brane configurations^, provided we are regarding 
one of the canonical parameters to be special. 

Notice also that the study of F-theory compactifications^ is effected by looking at elliptic 
fibrations described by elliptic curves of the Tate form similar to the one of Eq. fllGOp where 
now the coefficients are polynomials in the coordinates of the base manifold B. An important 
limit studied in these investigations is the Sen limit^ or orientifold limit. This is achieved by 
demanding that the axio-dilaton to be constant almost everywhere in type IIB space-time. 
Explicitly this limit is realized by setting 

in the Tate-form. In the two center stw-context aj are polynomial invariants of homogneous 
degree 2i and ag = 0. For the stu model we have = —LM, = D. For the st^ and 
models L = hence = D and a corresponding limit is achieved by sending the invariant 
D to zero by e.g. scaling down one of the charges to zero. This is precisely what happened 
in Section V.B. and also in the situation described by Eq. (11481) . 

Finally our triality symmetric Seiberg- Witten curve of Eq. (1186P also makes its presence 
in the F-theory contexlr^. Here the basic idea is to deform away from the special point in 
moduli space where the orientifold picture applies. This important deformation is effected 
by a curve of the form = + f{u)x + g{u) where / and g are polynomials in u of 
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degree two and three. Sen has shown that this curve can be cast in the form of Eq.f ll86p 
provided we make a suitable mapping of the parameters involved in the two different physical 
contexts. Notice that in the original papers the relevant curves^ii^ were not parametrized 
by the four-qubit invariants H, L, M, N and F as in Eq. (11861) . Instead in these papers four 
parameters were employed, which are easily identified with the four canonical parameters of 
four-qubit states of Eq. (!32|) . In particular in the F-theory context the relevant deformation 
is a one corresponding to splitting of the six coincident zeros of the discriminant away 
from each other, i.e in the orientifold picture moving the four coincident seven branes away 
from the orientifold plane. The collective coordinates of the background describing such 
a physical situation are that of an = 1 supersymmetric 5*0(8) gauge theory in eight 
dimensions with moduli space characterized by a complex scalar field $ . At a generic 
point in moduli space the vacuum expectation value of this field has the form similar to the 
canonical form of the central charge of Eq. (ll65p we used in Section V. in the different context 
of = 8 supergravity. Now one can check that the four canonical complex parameters 
ci,C2,C3 and C4 of ($) are behaving exactly like the canonical parameters of a four-qubit 
state. Again the polynomial expressions used by Seiberg and Witten in their curve of the 
Eq. (ll86p form are the ones for the four-qubit invariants H, L, M, N and F now with the 
simple identification's Cj = mj. Our four-qubit analysis as presented in this paper shows 
similar symmetry properties (permutation symmetry of the stu truncation connected to 
triality of 5*0(8) within A^ = 8 supergravity ). This might indicate that the physics of 
two-center black holes could be another arena where this curve plays a basic role. Of course 
these structural similarities could be superficial however, in any case the possible physical 
ramifications should be explored further. 

VII. CONCLUSIONS 

In this paper we have shown how the ^/-duality invariants introduced by Ferrara et.al.- 
characterizing two-center extremal black hole charge configurations in the stu, st^ and 
models oi N = 2, d = A supergravity can be understood as entanglement invariants of four- 
qubit systems. In this entanglement based picture the geometric and algebraic meaning of 
these invariants is displayed in a nice and unified manner. For one of the entanglement 
invariants that have not yet made its debut to the supergravity literature we have found a 
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distinguished role. It is the hyper determinant of type 2x2x2x2 which is the generahzation 
of Cayley's hyperdeterminant featuring the macroscopic black hole entropy formula in the 
stu model. For the special example of the BPS DODA — D2DQ split in the model we have 
demonstrated that this polynomial invariant of order 24 seems to govern important issues 
of consistency for the two center solutions. 

We have also introduced a quartic polynomial featuring the algebraically independent 
four-qubit invariants. For our simple example we have shown that the property that this 
polynomial has real roots provides a necessary condition for the consistency condition to 
hold. The resolvent cubic of this fundamental polynomial can be cast into a cubic of Weier- 
strass canonical form which in turn can be used to define an elliptic curve associated to 
the two-center charge configuration. After switching to the more convenient Tate form of 
this curve we have shown that this association is natural, meaning that the Tate form is 
displaying combinations of the algebraically independent invariants of physical meaning. 
Indeed, the Tate form falls short of terms step by step as we perform truncations to the 
st"^ and models, restrict attention to mutually local charge configurations, or perform the 
degenerate limit resulting in splits into small black holes. The discriminant of this elliptic 
curve is just the hyperdeterminant, a quantity also featuring the j invariant of the curve. 
For our example the structure of the j invariant nicely encapsulates the basic properties of 
the DODi - D2D6 split. 

The mapping from two center black hole charge configurations to elliptic curves is many 
to one. Moreover, for the stu model only a special class of curves with ae = in Eq.f ll30p 
can be reached. We observed that the vanishing of can be attributed to the vanishing 
of a polynomial of degree 12 characterizing the stu model. Based on this we conjectured 
that our N = 2 picture can be substantially generalized also incorporating two center charge 
configurations of = 8 supergravity. Here Cartan's quartic invariant and its polarizations 
should make their presence together with the usual symplectic invariant and a new invariant 
of order 6 introduced recently—. The discriminant of this curve of more general type then 
would play the role of a generalization of our hyperdeterminant, with probably similar 
physical meaning than its stu descendent. 

Finally we presented other physical contexts where our four-qubit invariants parametriz- 
ing elliptic curves also play an important role. Here we have given a new look to the triality 
invariant curve originally introduced by Seiberg and Witten^^. In its new form this curve 
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is featuring four-qubit invariants also displaying permutation invariance. In this new set- 
ting we have also invoked the F-theory interpretation of this curve as was given by Sen^^. 
Taken together with the two-center black hole context studied in this paper the main uni- 
fying theme in these scenarios seems to be the presence of a hidden torus and its extra 
set of PSL(2, Z) modular transformations. Note that in our investigations elliptic curves 
appeared merely as natural mathematical objects nicely encapsulating the information on 
the structure of two-center [/-duality invariants. However, apart from this we have also 
presented some evidence that the hyperdeterminant and the j invariant associated to these 
curves might contain useful information on issues of marginal stability, domain walls and 
split attractor flows. The physical background underlying the use of elliptic curves in the 
F-theory and N = 2 supersymmetric SU{2) gauge theory context is well-known. Is there 
any deeper physical reason also accounting for the occurrence of elliptic curves related to 
the structure of multi-center black hole solutions? 
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